Spin tunneling and topological selection rules for integer spins 
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as M1112 (Refs. |6j andO) and FegBEl The corresponding- 
classical anisotropy energy E Ztn (0,(f>) = — As 2 cos 2 9 + 
B n s n sin™ 9 cos(n0) has the shape of a double-well po- 
tential, with the easy axis pointing along the z di- 
rection. It is obvious that the anisotropy energy re- 
mains invariant under rotations around the z axis by 
multiples of the angle r\ = 2n/n, i.e., lZ ZtV E z (9, (f>) = 
E z (9, 4> + -q) = E z (9, (ft). It is convenient to describe the 
symmetry of the anisotropy energy by the cyclic group 
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We present topological interference effects for the tunneling of a single large spin, which are caused 
by the symmetry of a general class of magnetic anisotropies. The interference originates from spin 
Berry phases associated with different tunneling paths exposed to the same dynamics. Introducing a 
generalized path integral for coherent spin states, we evaluate transition amplitudes between ground 
as well as low-lying excited states. We show that these interference effects lead to topological 
selection rules and spin-parity effects for integer spins that agree with quantum selection rules and 
which thus provide a generalization of the Kramers degeneracy to integer spins. Our results apply 
to the molecular magnets Mni2 and Fes- 
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In recent years the tunneling of spin has attracted 
much attention since experiments provided strong ev- 
idence for the existence of spin tunnclinc-ia molecu- 
lar magnets, such as Mma-acetate (Mn^lu'tHa and Fes- 
triazacyclononane (Fes).aLi In Refs. [j], ^], [| |] the relax- 
ation time t(H z ) of the magnetization as a function of 
the longitudinal magnetic field H z was measured, with 
t{H z ) exhibiting large peaks due to thermal-assisted tunp. 
neling of the spinBEl In Rcf. g, the incoherent Zener ratecl 
as a function of the transverse magnetic field H x was 
measured, from which quantum oscillations of the tun- 
nel splitting E mm i(H x ) can be ohserved.El In a coherent 
spin-state path-integral approach,tl these quantum oscil- 
lations and their associated spin-parity effects in the case 
of quadratic anisotropies can be viewed as a result of in- 
terfering Berry phases carried by spin tunneling paths of 
opposite windings, which arc modified in the pres- 
ence of a transverse-field H x (Ref. [l2|) or crystal fields of 
cubic symmetry with multiple degenerate minima. t3 

The results obtained in the present work can be sum- 
marized as follows. If the Hamiltonian for a single spin 
s> 1 contains a magnetic anisotropy of arbitrary order 
n that is transverse to the easy axis and symmetric un- 
der rotations of 27r/n around the easy axis, there are n 
different paths that possess the same dynamical action 
but different spin Berry phases and which thus interfere 
constructively or destructively. Using a spin path inte- 
gral that we generalize to excited states we will show that 
the tunneling of the spin between ground or low-lying ex- 
cited states that are doubly degenerate, \m) and |— m), 
|m| < s, is completely suppressed if n is not a divisor of 
2m (i.e., 2m/ n ^ Z); this result represents a generaliza- 
tion of the Kramers degeneracy to integer spins s. We 
also provide the connection to quantum selection rules 
using standard spin algebra, and we discuss the effect of 
magnetic fields on these selection rules. 

We consider a general single-spin Hamiltonian TL z _ n = 



-l)i7>fte»,n»7 = where 



T^a,ip performs rotations about the a-axis, a — x,y, z, by 
the angle ip, and X denotes the identity operator. 

For the following calculation it proves favorable to 
choose the easy axis along the y direction. Then the 
spin Hamiltonian and the anisotropy energy are changed 
into 
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-AS 2 v + B n {S n + + S n _), 



where now S± = S z ±iS x , and 
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y,n\ 



-As 2 sin 2 9 sin 2 <f> + B n s n [(cos 9 
-i sin 9 cos </>)" + (cos 9 — i sin 9 cos 



(1) 



<t>n (2) 



We note again that C y ^ n E y ^ n = E Vtn , which will play a 
central role in what follows. 

We are interested in the tunneling between the eigen- 
states |m) and \—m) of —AS 2 , corresponding to the 
global minimum points (9 = n/2,<j) = —ir/2) and (9 = 
it/2, cj) = +7r/2) of E y , n . For this we evaluate the imagi- 
nary time transition amplitude between these points. For 
the ground-state tunneling (m — s ^> 1) this can be-done 
by means of the coherent spin-state path integralHEJ 



7r/2,+7r/2 



vn e~ SE , 



(3) 



7r/2,-7r/2 



-AS 2 + B n (Sl + S 7 1), with easy-axis (A) and trans- sure of the S 2 sphere, and S E = /<f dr[is^(l 



where (3 — 1/fcgT is the inverse temperature, T>fl = 
Ii T dVL T , dVl T = [Air /{2s + l)]d(cos9 T )d<p T the Haar mea- 





verse (B n ) anisotropy constants satisfying A B n > 
0. Here, n is an even integer, i.e., n — 2,4,6,..., 
so that Ti. z is invariant under time reversal. Such 
Hamiltonians are relevant for molecular magnets such 
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Ey t , n ] the Euclidean action, where the first term in. 
Se defines the Wess-Zumino (or Berry phase) term,Efl 
which gives risp to topological interference effects for spin 
turmelrngBBBB 
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FIG. 1. Interference between n = 6 paths (small arrows) 
moving from — y to y direction (large arrows), which is de- 
structive if n is not a divisor of 2m. The grey surface area 
<3?(m)/m depicts the Berry phase $(m) of the upper right 
path. Note that there is no simple pairwise interference; in- 
stead, all six paths interfere with each other [see Eq. (g)]. 

For the tunneling between low-lying excited states 
(1 <C \m\ < s) we now derive a generalized coherent 
spin-state path integral. It is knownE-3 that one can con- 
struct a coherent state system {T \ipo)} out of any start- 
ing vector \ifjo) of a Hilbert space H, with T being a rep- 
resentation of any Lie group operating on this EL Thus, 
applying this to our case, we can choose any spin state 
\m) el s = {\m) I — s < m < s} as starting vector to- 
gether with the Lie group SU2 , leading to the generalized 
coherent spin states 



\S7,m) 



-iS y 9 e -iSzX 



(4) 



where <fi, 8, and \ are the Euler angles. Note that the 
dispersion of S 2 is not minimal for \m) 7^ \s). How- 
ever, as long as 1 \m\ < s, the spin fluctuations re- 
main small. By means of Schur's Lemma, one can prove 
as in Ref. [l5] that the integral J dO \fl,m) (f2,m\ oc 
T. So we have to calculate one single diagonal ele- 
ment in order to fix the proportionality constant. For 
this it turns out to be convenient to express the states 
1 17, m) by 2s spin-1/2 states that are individually ro- 
tated, i.e., |f) = cos(0/2) |T) + e** sin(0/2) |j>, = 



e-^sin(6»/2) |t) - cos(0/2) ||), where the sing 
states are represented in the north pole gauge, 
we obtain 



i valued 
3 Thus, 



\n,m) 



2s 



-1/2 



{jk} 



ii<—<J s - 
i<ife<2» 



where \ denotes the s — m rotated down spins \ij k ), 
at the positions jk ■ Evaluation of the amplitude 

(f2,m\s) = (n,m\ hh---hs) 




s+m 



sin- 



ce) 



leads to / dfl {s\f2, m) {f2, m\s) = 4?r/(2s + 1). Hence, 
the integration measure for all \f2, m) states is the same. 
Next, to derive the path integral, we need the quan- 
tity (f2,m\S\ n,mj = (mfl + Oyi)) {fi,m\n,rnj^ 

where we must choose the value of m close to s. Con- 
sequently, we have to evaluate the overlap between two 
states separated by the angles 89 = 9 — 6 and 8(f) = <j>—<j>, 

///|J?. in ) = (terms containing (f' | J.') and (j' | f')) 
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2 2 



sin - sin -e 
2 2 
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First-order expansion in 89 and 8(f) leads to the 
overlap between two infinitesimally separated states, 

(^S7,m\f2,mJ w 1 + im(cos9 — l)8<f>. After slic- 
ing the imaginary time, .interval f3 into N pieces 
of length e = /3/JV,Edii9 the transition amplitude 
(12,m(r„ +1 ) |1 — eH\ f2,m(r n )) between two time steps 
T n+ i and T n = ne can be approximated by (1 — 
eTL[mfl\) {fl,m{T n+ i)\fl,m(r n )), where H[mf2] is the 
diagonal element of the Hamiltonian. Collecting the var- 
ious terms, we obtain 

f2,m\e-^ H \ (2,m}= f VVL e^C"*)-^, ( 7 ) 



where <Sd yn = f At 7i.[mf2] is the dynamical action 



and <E>(m) = m J^ 2 



/2' 



(1 — cos0) the Berry phase, with 



<&(m)/m measuring the area of the spherical triangle be- 
tween the tunneling path and the great circle through 
the north pole of S 2 (see Fig. 0). 

Now we possess the tool to compute the path integral 
in Eq. (|J) also for low- lying excited states |m), m < s. 
Since substitution of s by m in Eq. (|^) does not change 
the symmetry of E v>n , Sdyn remains invariant under C v ^ n 
operations, whereas lZ y . v &(m)/m = $(m)/m + A-k/ti. 
Thus,we can divide the S 2 surface area of integration in 
Eq. (0) into n equally shaped subareas A n , so that we 
can factor out the following sum over Berry phase terms 
without the need of evaluating the dynamical part of the 
path integral (i.e., the following result is valid for all m), 



oc yv( 2 *v«)( 2 fc-i)m = 

k=l 



(5) 



sin(27rm) 
sin[(27r/n)7r 

(8) 
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which vanishes whenever n is not a divisor of 2m. Thus, 
the tunnel splitting energy E m .- m between the states |m) 
and | — m) vanishes if 2m/ 'n ^ Z. However, if 2m/ n G Z, 
the variable m must be extended to real numbers, i.e., 
m — > fi € R, in order to calculate the limit zx — ► m of 
the ratio of the sine functions, which is plotted in Fig. || 
for a special case. In order to visualize the interference 
between the Berry phases in Eq. (||), we select one rep- 
resentative path of each subarea A n . Then the vanishin 
of the amplitude in Eq. (g) for the case n = 6 (see Fig 
can be thought of as a destructive interference between 
six different paths. Note that there is no simple pairwise 
cancellation .E!l Also, it is important to note that from the 
semiclassical point of view the total classical energy of the 
spin system E must be conserved during the tunneling 
process. Therefore the semiclassical paths do not follow 
the local minima of E with respect to 9 alone, except-fpt 
the case of purely quadratic (i.e., n = 2) anisotropies.t3't2l 



sin (ll2n) /sin (^2jr/n) 




FIG. 2. Sum over Berry phase terms for s = 10 and n = 4. 

From the above treatment we conclude that tunneling 
between two degenerate spin states |m) and |m') = \—m) 
is topologically suppressed (i.e., the twofold degeneracy 
is not lifted) whenever n is not a divisor of 2m. Since n 
is even this excludes immediately tunneling for all half- 
odd integer spins s (for all m and n) , in accordance with 
Kramers degeneracy. For the s integer, however, tunnel- 
ing can be either allowed or suppressed, depending on 
the ratio 2m/ n. In the latter case, the twofold degener- 
acy of spin states is not lifted by the anisotropy, and we 
can view this result as a generalization of the Kramers 
theorem to integer spins. For the special case n = 2, we 
immediately recover the spin-parity effect found before.EHI 

Now let us take a look at the tunnel splitting energy 
E m m' by means of higher-order degenerate perturbation 
theory in terms of a standard spin operator formalism. 
Applying the resolvent techniqueQ we find 



(m— m') / n 
tin 



Yll™-™'^- 1 A [m? - {m- knf] 



^m — m'yjn j 



(9) 



where j E Z. Note that E m m' vanishes if 2m /n Z 



7~ty,n divides the Hilbert space H s into a direct sum of 
invariant subspaces hk = {|m') | s — m' = k modrt}, i.e., 
M s = ©fc = i hk- Thus, the topological interference effects 
due to the Berry phase presented above correspond to 
these quantum-mechanical selection rules. 

Next we take the external transverse magnetic field H z 
for the case n = 2 into account. For this, we consider the 
following Hamiltonian {A! =A- 2B 2 , B' = 4B 2 ): 



H 



y,2 



-AS* + B'S; + h z S z + Am z C, (10) 



yielding the classical energy E y ^{u) = 
—A'm 2 [(l — w 2 ) sin 2 (f> — C] + B'm 2 u 2 + h z mu, where 
h z = gfisHz and u — cosO. The constant C — 1 is de- 
termined by the boundary conditioncl u(<j> = ±7r/2) = 
for h z = 0. By completing the square in E y ^(u) by the 
substitution u = u + h z /2m(A' sin 2 <f> + B'), the Jacobian 
of which is unity, we end up with the anisotropy energy 



E y , 2 {u) = -A'm 2 [(1 - u 2 ) sin 2 <f) — i\ + B 
1 h 2 



'„,2 ~1 

m ,u 



4A> 



B> 



(11) 



E y ^(u) and the new integration limits for u are invariant 
under C y ^ operations. Therefore, we can factor out two 
exponential terms with the following Berry phases: 



<i>±(m) — ml 

J — 7T 



7r/2,-37r/2 

( 

/2 



1 



2m(A' sin 2 + S') 



±7rm 1 



2my/B'(A' + B') 



(12) 



for 



m = — m. 



This is a consequence of the fact that 



from which we directly obtain (— 7r/2 |e P n y- 2 1 + tt/2) cx 
cos($ + ) and thus, Berry phase oscillations in H z with 
period 2y/B'(A' + B') = 2y/4B 2 (A + 2B 2 ), and with 
phase shifts \Q±(m) — $>±(m— 1)| = n, regardless between 
which states |m) and \—m) the tunneling takes place. 
Since < 1$+ — $_| < Avm, there are 2m zeroes of the 
tunnel splitting £ , m ._ m .E3 These findings agree with the 
exact diagonalization of Eq. ([To]) and have been observed 
in Fe8.El Note that these interference effects are due to 
the symmetry C v ^ 2 . Hence, when breaking this symme- 
try with a magnetic field that is tilted away from the z 
axis (in the xz plane) by the angle ih- the oscillations be- 
come suppressed with increasing i/jfB in agreement with 
Rcf. ||. In contrast to previous work on n = 2 (Rcf. |T^ ) 
we did not use approximations based on semiclassical 
paths (given by energy conservationXpr on semiclassical 
Wentzel-Kramers-Brillouin methodsEj to find the zeroes 
of the tunnel splitting. 

Next, we consider the case n = 4, 

H yA = -AS 2 + Bi {S\ + Si) + h z S z + Am 2 C, (13) 

and its corresponding classical energy 

E yA = -Am 2 [(1 - u 2 ) sin 2 <j>-C]+ 2B 4 m 4 [u 4 - 6u 2 

x (1 - u 2 ) cos 2 (/)+(!- u 2 ) 2 cos 4 <j>] + h z mu. (14) 
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Again, the boundary conditions determine the constant 
C = 1 for h z — 0. Unfortunately, it is no longer possi- 
ble to complete the quartic polynomial in Eq. ([Ti]). The 
same problem persists for n > 4. This means that in 
contrast to n = 2, the C y ^ symmetry of the anisotropy 
energy E y ^ without transverse field H z cannot be re- 
stored anymore. However, there still remains the reflec- 
tion symmetry between the conjugated paths 7+ = (9, <p) 
and 7_ = (8, —(f>) that wind around the z axis counter- 
clockwise and clockwise, respectively. The topological 
part of the main contribution to the tunneling ampli- 
tude exhibits oscillations, as shown in Fig. [| This part 
comes from the two semiclassical paths u c i,±(0) that sat- 
isfy E Vt 4 — and whose Berry phases, defined by 

/7r/2,-3Tr/2 
d<P{l - Re{u cl>± }) , (15) 
-7T/2 

interfere with each other. Thus, the level splitting in- 
duced by tunneling shows quantum oscillations as a func- 
tion of H z due to this oscillating Berry phaseEj Again, 
for a magnetic field with ip ^ kn/2, k <G Z, the reflec- 
tion symmetry between conjugated paths is broken (see 
above), and therefore the more if> is detuned from kiri2. 
the more the Berry phase oscillations are suppressed.tJ 
Note that the level splitting contains m zeroes, m = 
8,9,10, in the interval -4 T< H z < 4 T (see Fig. |), 
and that the Berry phase oscillations are shifted by tt/2 
when the parity of m changes (see m = 8, 9, lAin Fig. 0), 
all in agreement with exact diagonalizationEd Thus, we 
have obtained a new spin-parity effect for n — 4 in the 
case of integer spins s which results from Berry phase os- 
cillations. This effect should be observable for tjj = kir/2 
in molecular magnets such as Mni2- Since for all n two 
conjugated paths (see above) can be found, Berry phase 
oscillations and their associated parity effect are max- 
imal whenever ip = fc27r/n.E3 This provides an exper- 
imental means to determine the order n of transverse 
anisotropies. 




H Z [T] 



FIG. 3. Topological part of the tunneling amplitude for 
n — 4 and s = 10 as a function of H z [see Eq. Jlq)]. Solid 
line: m = 10. Dashed line: m = 9. Dotted line: m = 8. The 
right boundary of cos($ c i) is determined by < 2|$ c i| < 47rm. 
The values A = 0.54 K and Ba = 8.5 x 10 -5 K approximate 
roughly the anisotropies in Mni2 (Ref. 



In conclusion, we have shown that the tunnel splitting 
vanishes due to interfering spin Berry phase and thus 
the level degeneracy is not lifted whenever n is not a 
divisor of 2m. This generalization of the Kramers degen- 
eracy to integer spins holds not only for the ground state 
but also for low-lying excited states m and for transverse 
anisotropies of arbitrary order n; the presence of mag- 
netic fields leads to oscillations of the tunnel splitting. 
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